Ballistic conductance of Ni nanowire with a magnetization reversal 
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The approach proposed by Choi and Ihm for calculating the ballistic conductance of open quantum 
systems is generalized to deal with magnetic transition metals. The method has been implemented 
with ultrasoft pseudopotentials and plane wave basis set in a DFT-LSDA ab-initio scheme. We 
present the quantum-mechanical conductance calculations for monatomic Ni nanowire with a single 
spin reversal. We find that a spin reversal blocks the conductance of d electrons at the Fermi energy 
of the Ni nanowire. On the other hand, two s electrons (one per each spin) are perfectly transmitted 
in the whole energy window giving 2Go for the total conductance. The relevance of these results in 
connection with recent experimental data is discussed. 



Understanding electron conduction through atomic 
and molecular wires connecting two macroscopic elec- 
trodes is presently a very active research topic. Con- 
ductance measurements for metal nanocontacts display 
flat plateaus and abrupt drops during the elongation at- 
tributed to atomic rearrangements. The last conductance 
step before breaking, corresponding to a monatomic 
nanowire, has a value 2Gq (Go = e 2 /h is the conductance 
quantum per spin) in monovalent metals such as gold 
which can be rationalized with the free propagation of 2 
(1 per each spin) valence s-electrons. Surprisingly, a low 
value of the last conductance step was also observed in 
nanowires and nanocontacts of transition metals having 
partially occupied s and d states. Oshima et al. @, who 
worked in vacuum, at variable temperature, and with the 
possibility of a magnetic field, found a minimal conduc- 
tance step in Ni nanocontact preferentially near 2 and 
4 at RT and zero field, near 4 at 770 K and zero field, 
and near 3 (occasionally near 1) at RT with an applied 
magnetic field. Ono et al. [3j, reported again 2 for Ni 
in zero field, and 1 for Ni in a field. Recently, Ugarte et 
al. Q observed 1 conductance quantum at room tem- 
perature and zero magnetic field in atomic chains (the 
ultimate configuration before the nanocontact rupture) 
of transition metals such as Co and Ni. 

Several theoretical methods are available to study 
the transport properties of atomic scale conductors. 
Among these, there are methods based on nonequilib- 
rium Green's functions combined with a localized basis 
set |H IS- A method based on a system-independent 
wavelet basis set to calculate the Green's function 
has been also presented recently. Another approach for 
conductance calculation relies on the direct solution of 
the scattering problem for an open system and on the 
application of the Landauer-Buttiker formula relating 
the ballistic conductance to the total transmission coef- 
ficient at the Fermi energy. The scattering approach has 
been applied by Lang and co-workers Q, and Tsukada 
and co-workers using the jellium model for the elec- 
trodes. The layer KKR approach has been used to study 



the spin-dependent tunneling in magnetic tunnel junc- 
tions Lately Wortmann and co-workers formu- 
lated a transfer matrix approach based on computation 
of single-electron Green function in the linearized aug- 
mented plane wave (LAPW) basis. 

We are interested in a formulation based on plane 
waves and ultrasoft pseudopotentials, a calculational 
technology which we are applying to other problems in- 
volving transition metals j 1 21 ] - Recently, Choi and Ihm 
[liij l proposed a first-principles method based on plane 
waves to solve the scattering problem with real atomic 
contacts. In their formulation atoms are described by the 
Kleinman-Bylander-type pseudopotentials. We general- 
ized this scheme to deal with ultrasoft pseudopotentials 
[lij in order to apply it to magnetic transition metals. 

In a recent paper |l5| we showed that a magnetization 
reversal (sharp domain wall) built inside a Ni nanowire 
effectively blocks the d channels allowing only two s elec- 
trons (one per each spin) to propagate and thus leaving 
only two conductance channels. That result was argued 
for an infinite monatomic Ni wire with a periodically 
repeated spin reversal by simply counting the number 
of bands crossing the Fermi level. However, no actual 
quantum-mechanical calculation of conductance was car- 
ried out. Now we substantiate our earlier suggestion by 
calculating the ballistic conductance of monatomic Ni 
nanowire with a single spin reversal. Wc present here 
the intrinsic conductance of monatomic nanowires due 
to only the scattering on the magnetization reversal and 
do not consider any additional contact resistance arising 
from the tip-wire junctions. 

METHOD 

We study electron transport in the linear regime (in 
vanishingly small applied bias) in the open quantum sys- 
tem consisting of a scattering region (0 < z < L) at- 
tached to the left (z < 0) and to the right (z > L) 
to semi-infinite generic electrodes. We assume that the 
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electrons move ballistically in the self-consistent poten- 
tial and all elastic scattering takes place in the scattering 
region which can be a monatomic wire, a molecule or just 
a single atom. The self-consistent potential is obtained 
performing ground state DFT calculations with a super- 
cell containing the scattering region and sufficiently large 
pieces of left and right electrodes to get the bulk-like po- 
tential for the electrodes. A wave propagating at energy 
E with in the xy plane incident from the left electrode 
on the scattering region will form a scattering state f ki 
which due to the supercell geometry has a Bloch form in 
the xy direction: 

* k± (r ± +R ± ,z;a) = e lk ^4v(r i ,z;a), (1) 

where R± are wire supercell lattice vectors. The spin- 
polarized systems are treated within the local spin- 
density approximation (LSDA) . Therefore electrons with 
different spin directions move independently in different 
self-consistent potentials and the scattering problem can 
be solved separately for the two spin polarizations. In 
the same way, different k^ do not mix and can be con- 
sidered separately. We will therefore omit both k^ and 
a from now on. Note that if the supercell is sufficiently 
large in the xy plane one can limit the calculation to the 
two-dimensional T point, k^ = 0. 

We describe the atoms of transition metals such as Ni 
by ultrasoft pseudopotentials Q . In the ultrasoft pseu- 
dopotential scheme the stationary state 'J at the energy 
£ is a solution of single-particle Kohn-Sham equation 
which as in Ref.|lfil| we write in the form: 

[-V 2 + Kff(r)] *(r) + ^ D^M^Piir-Rj) = £?*(r). 

Iran 

(2) 

with Din = DLn-Eq'mn- Th e coefficients D' mn and ^ mn 
characterize the pseudopotential which is constructed us- 
ing the set of projector functions (3^ associated with 
atom I (see Refs.[lJ,[l3). 

Deep within the electrodes the scattering state orig- 
inating from the rightward propagating Bloch wave tpj in 
the left electrode has an asymptotic form: 

!ipj + J2 njipi, z < 
ieL 
2J Ujipi, z > L 

where summation over i € L (i G R) includes the general- 
ized Bloch states in the left (right) electrode at energy E 
which propagate or decay to the left (right). The general- 
ized Bloch states (both propagating and decaying) con- 
stitute the so-called complex band structure of a solid. 
In Ref. we presented in detail the calculation of the 
complex bands with ultrasoft pseudopotentials. Having 
obtained the complex bands for the electrodes we calcu- 
late the transmission and reflection coefficients tij and 
by generalizing in a straightforward way the approach by 



Choi and Ihm [13| to the case of ultrasoft pseudopoten- 
tials (the details will be presented elsewhere). 

The ballistic conductance G is related to the total 
transmission T at the Fermi energy by the Landauer- 
Biittiker formula G = GqT (Go = e 2 /h is the conduc- 
tance quantum per spin). The total transmission is given 
by: 

T = ^|7y| 2 = Tr[T+T], (3) 

ij 

where T is the matrix of normalized transmission ampli- 
tudes Tij — y/Ii/Ij ■ tij and Ij is the probability current 
of the state tpj in the z direction. Note that only right- 
ward propagating states in both left and right electrodes 
should be considered so the matrix T is of dimensions 
MrMl where Ml and Mr are the number of propagat- 
ing modes in the left and right electrodes, respectively. 

The eigenvectors of the Hermitian matrix T + T deter- 
mine the coefficients of a unitary transformation from 
the set of Bloch states tpj to the nonmixing conductance 
eigenchannels [l^. I n the eigenchannel basis the matrix 
T + T is diagonal. Calling its eigenvalues TJ, the conduc- 
tance is a sum of independent contributions from each 
eigenchannel: 

G = G J2Ti, (4) 

i 

where Ti gives the transmission probability for z-th eigen- 
channel. 

In order to calculate the total transmission T one 
needs to know the current Ij carried by the propagat- 
ing Bloch state tpj in the z direction. The tpj is actu- 
ally pseudo-wave function which coincides with the all- 
electron wave function only outside the core regions. The 
true current does not coincide with the pseudo-current, 
and their difference is not trivial. We recently showed 
how the pseudo-current should be augmented to obtain 
the real current [l9j and that method is implemented in 
the present calculation. 

The calculation of ballistic conductance then proceeds 
as follows. First, we perform the supercell DFT elec- 
tronic structure calculation with the plane-wave code 
(PWscf ) [i(J to get the self-consistent potential V c g and 
the screened coefficients D^, For spin-polarized sys- 
tems both V e ff and D I mn will depend on the spin direc- 
tion. Second, we calculate the complex band structures 
of the left and right electrodes. The unit cells of the 
electrodes are chosen within the supercell but far enough 
from the scattering region where the potential becomes 
bulk-like (for details concerning complex band structure 
calculations see Ref. ^(| ) . At last we calculate as in Choi 
and Ihm's work the transmission coefficients £y for each 
rightward propagating state tpj in the left electrode and 
obtain the total transmission T and eigenchannel proba- 
bilities Ti by diagonalizing T + T as described above. 
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RESULTS 

We apply the scheme outlined in the previous section 
to calculate the ballistic conductance of monatomic Ni 
nanowire with a single spin reversal. No contact scatter- 
ing resistance is considered here so that the left and right 
electrodes are merely monatomic wires magnetized in op- 
posite directions. We consider a wire the left (right) half 
of which has a positive (negative) magnetization. A su- 
percell containing 12 atoms (6 with up and 6 with down 
magnetization) was used to calculate the self-consistent 
potential and the total magnetization was constrained 
to be zero (for computational details see Ref . |lfij ) . The 
scattering region and the unit cell of the left and right 
part of the wire are shown schematically in Fig.l. Sym- 
metry requires that the total transmission be the same 
for electrons of both spin directions. Therefore we chose 
to consider only electrons with down spin polarization 
incident from the left on the spin reversal region. The 
complex bands of the left (right) electrodes needed for 
the transmission calculation are those shown on the right 
(left) panel in Fig. 2 of Ref. [3. 

In Fig. la and lb we present the total transmission 
and the contributions from each band, respectively, as a 
function of energy for the Ni nanowire with a spin re- 
versal. A propagating state in the left part of wire can 
propagate through the spin reversal and contribute to the 
total transmission only if it can be matched with the ap- 
propriate propagating state to the right. For example, at 
the Fermi energy there are 6 propagating states in the left 
part of the wire. However, only two s — d z 2 states have 
nonzero transmission (Fig. lb) since in the right piece 
of the wire there is only one s — d z 2 propagating state. 
The d states can be matched only with decaying states 
in the right part of the wire and are therefore completely 
reflected by the spin reversal. At an energy of about 
E = Ep — 0.6 eV the two d xz , yz states start contributing 
to the total transmission because there are two states of 
the same symmetry in the right part of wire. Note that 
the d XV:X 2_ y 2 states can not get through the spin reversal 
at any energy and never contribute to conductance. 

Fig. lc shows eigenchannel transmissions T,. The two 
s — d z 2 states are combined to form a single perfectly 
propagating conductance channel (mostly s like). An- 
other channel (mostly d z 2 ) orthogonal to the first appears 
at energy about 0.6 eV below the Fermi energy. It is quite 
narrow and has a smaller transmission probability (about 
0.9 at the maximum). The two d xz ,yz states form two 
eigenchannels of the same symmetry. They contribute to 
the total transmission between 0.6 eV and 2.2 eV below 
the Fermi energy and have almost unit transmission at 
the maximum. 

The total ballistic conductance for a wire subject to 
left and right chemical potentials /i^ and fin (/i^ > /ir) 
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FIG. 1: The electron propagation through a spin reversal 
in a monatomic Ni nanowire: a) total transmission; b) con- 
tributions from each band; c) transmission eigenvalues T%. 
The bands and conducting channels are labeled by their main 
atomic character. The Fermi energy is chosen to be zero 
(dashed vertical lines). 



is: 

G = G ° V r L T a (E)dE. (5) 

ml - m ^ J m 

For hl — ► fJ-R, when G = Go'^2 (T T <t (Ef), we obtain 
an ultimate conductance of 2Go for a monatomic Ni 
nanowire that disagrees with the experimental observa- 
tion of a value rather closer to Go 0. We should stress 
however that the exact equivalence of s channels for the 
two spin polarization, leading to the factor 2 in our re- 
sult, is likely to be an artifact of our chosen symmetrical 
geometry. In reality, there will generally be no up-down 
spin symmetry, and there might also be a more complex 
magnetic configuration. It seems possible that in such 
more complex geometries one of the s spin channels could 
be blocked as well, which would lead to Go as observed. 
We plan to elaborate further on this line in later work. 
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CONCLUSIONS 

We generalize to magnetic transition metals the ap- 
proach of Choi and Ihm for calculating the ballistic con- 
ductance of an open quantum system. The method has 
been implemented with ultrasoft pseudopotentials and 
plane wave basis set in a DFT-LSDA ab-initio scheme. 
We applied our method to calculate the ballistic conduc- 
tance of monatomic Ni wire with a single spin reversal. 
We found that at the Fermi energy there is only one con- 
ductance channel per spin (s like) formed by two s — d z 2 
propagating states while all d electrons are completely 
reflected by a spin reversal as previously guessed from 
band structure calculations ^5|. The s-like channel has 
in fact a unit transmission probability at any energy. 
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